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We probe the scale dependence of the primordial spectrum in the light of the three-year WMAP 
(WMAP3) alone and WMAP3 in combination with the other cosmological observations such as 
galaxy clustering and Type la Supernova (SNIa). We pay particular attention to the combination 
with the Lyman a (Lya) forest. Different from the first-year WMAP (WMAP1), WMAP3's pref- 
erence on the running of the scalar spectral index on the large scales is now fairly independent of 
the low CMB multipoles I. A combination with the galaxy power spectrum from the Sloan Digital 
Sky Survey (SDSS) prefers a negative running to larger than 2a, regardless the presence of low £ 
CMB (2 < £ < 23) or not. On the other hand if we focus on the Power Law ACDM cosmology 
with only six parameters (matter density Q m h 2 , baryon density Qbh 2 , Hubble Constant Ho, optical 
depth r, the spectral index, n s , and the amplitude, A s , of the scalar perturbation spectrum) when 
we drop the low £ CMB contributions WMAP3 is consistent with the Harrison-Zel'dovich-Peebles 
scale-invariant spectrum (n s — 1 and no tensor contributions) at ~ la. When assuming a simple 
power law primordial spectral index or a constant running, in case one drops the low £ contributions 
(2 < I < 23) WMAP3 is consistent with the other observations better, such as the inferred value 
of as- We also find, using a spectral shape with a minimal extension of the running spectral index 
model, LUQAS+ CROFT Lya and SDSS Lya exhibit somewhat different preference on the spectral 
shape. 
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I. INTRODUCTION 

It is commonly believed that the early Universe has 
experienced a stage of accelerated expansion which is 
known as inflation 0, 0, Q . During inflation one of the 
tiny patches of the very early Universe was supcrlumi- 
nally stretched to become our observable Universe to- 
day. The theory of inflation can naturally explain why 
the universe is flat, homogeneous and isotropic. Inflation 
is driven by the potential energy of a scalar field called 
infiaton and its quantum fluctuations turn out to be the 
primordial density fluctuations which seed the observed 
large-scale structures (LSS) and anisotropy of the Cos- 
mic Microwave Background (CMB) Radiation. Inflation 
typically ends with a period of reheating when the infia- 
ton starts rapid oscillations and it decays into standard 
particles. Then the universe starts the standard hot Big 
Bang evolution. 

Over the past decade, the theory of inflation has 
successfully passed several nontrivial tests. In partic- 
ular, the first year Wilkinson Microwave Anisotropy 
Probe(WMAP) llHESHIi has detected a large-angle 
anti-correlation in the temperature-polarization cross- 
power spectrum, which is the signature of adiabatic su- 
per horizon fluctuations at the time of decoupling ja] . 

The recently released three year WMAP (WMAP 3) 
[ToL ITTl Il2l ll3L Il4| has marked another milestone on 
the precision cosmology of CMB Radiation. The sim- 
plest six-parameter power-law ACDM cosmology, namely 



with matter density tl m h , baryon density flbh 2 , Hubble 
Constant Hq, optical depth r, the spectral index and 
amplitude of the scalar perturbation spectrum, n s and 
A s , turn out to be in good agreement with WMAP3 to- 
gether with the small-scale CMB observations such as 
BOOMERANG 0, the Arcminute Cosmology Bolome- 
ter Array Receiver (ACBAR [T^l. the Cosmic Back- 
ground Imager (CBI [l7|) and the Very Small Array (VS A 
111), LSS as measured by Two degree Field (2dF) 
and the Sloan Digital Sky Survey (SDSS) [2(J and with 
the Type la Supernova (SNIa) as measured by the Riess 
"gold" sample 12 ll and the first year Supernova Legacy 
Survey (SNLS) |22{. This agreement between the above 
"canonical" cosmological model and observations can be 
used to test a number of possible new physics, such as the 
equation of state of dark energy, neutrino masses, cosmic 
CPT violation, etc. 

Although the temperature-temperature correlation 
power (TT) of WMAP3 is now cosmic-variance limited 
up to £ ~ 400 and the third peak is now detected, the 
features discovered by the first year WMAP [|| Q are 
still present: the low TT quadrupole and localized os- 
cillating features on TT for I ~ 20 - 40 12]. Although 
the signatures of glitches on the first peak discovered by 
the first year WMAP have now become weaker, they do 
exceed the limit of cosmic variance |l2j . In fact, even be- 
fore the release of the first year WMAP Ref. [23| claimed 
oscillating primordial spectrum could lead to oscillations 
around the first peak of CMB TT power. These data 
have hence revealed many interesting detailed features 
which cannot be explained in the simplest version of the 
inflation model, and we are making much effort to con- 
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struct a realistic particle physics model of inflation. 

As for the shape of the primordial power spectrum, 
it is also noteworthy that a significant deviation has 
been observed by WMAP3 from the simplest Harrison- 
Zel'dovich-Peebles scale-invariant spectrum (n s — 1 and 
no tensor contributions), and that this feature is more 
eminent with the combination of all the currently avail- 
able CMB, LSS and SNIa (dubbed the case of "All" in 
|14|). On the other hand, a nontrivial negative run- 
ning of the scalar spectral index a s , whose existence 
was studied even before WMAP epoch .24] (for relevant 
study see also]25jL was favored by the first-year WMAP 
papers 0, 0, Q- " But its preference was somehow di- 
minished as corrections to the likelihood functions were 
made [H S HI] . Dramatically, the new WMAP3 data 
prefers again a negative running in the "All" combina- 
tion |14|. If confirmed, a nonvanishing running of a s 
would not only constrain inflationary cosmology signifi- 
cantly miMEJIilSHES&but also affect 
the cosmological constraint on the neutrino mass [3^ . l39| . 

In the present paper, using the Markov Chain Monte 
Carlo (MCMC) method, we aim to probe the scale depen- 
dence of the primordial spectrum in the light of WMAP3 
alone and that in combination with the other cosmolog- 
ical observations from galaxy clustering and SNIa. We 
pay particular attention to the combination with the Ly- 
man a (Lya) forest. 

We will show that, different from the first-year WMAP 
(WMAPI), WMAP3's preference on the running of the 
scalar spectral index on the large scales is now fairly inde- 
pendent of the low CMB multipoles £. We also find that a 
combination of WMAP3 with the galaxy power spectrum 
from SDSS prefers a negative running with more than 2a, 
regardless of the presence of low £ CMB (2 < I < 23). 
On the other hand, if we focus on the Power Law ACDM 
cosmology with only six parameters, we find that when 
one drops the low £ CMB contributions WMAP3 would 
be consistent with the Harrison-Zel'dovich-Peebles scale- 
invariant spectrum at ~ la. Assuming a simple power- 
law primordial spectral index or a constant running, in 
cases one drops the low £ contributions (2 < £ < 23) 
WMAP3 is consistent with the other observations bet- 
ter, such as the value of a&. 

On the other hand, using a spectral shape with a min- 
imal extension of the running spectral index model, we 
find LUQAS+ CROFT Lya gjj and SDSS Lya 
exhibit somewhat different preference on the spectral 
shape, although it is difficult to draw any definite con- 
clusion because the result depends on how we model the 
spectral shape. 

The rest of the paper is structured as follows. In Sec- 
tion II we describe the method and the data. In Section 
III we analyze the effects of low CMB multipoles on the 
determination of cosmolo gical p arameters using WMAP3 
DUES El 13 13, SNIa jHEl, 2dF and Sloan Dig- 
ital Sky Survey 3-D power spectrum (SDSS P(k)) [H, 
LUQAS+CROFT sample of Lyman a gjj and SDSS Ly- 
man a sample 4l| by global fittings using the MCMC 



technique. In Section IV scale dependence of fluctua- 
tions is analyzed with priors on the Hubble parameter, 
flbh 2 , and the cosmic age. 

Discussions and conclusions are presented in the last 
section. 



II. METHOD AND DATA 

To break the possible degeneracies among the varia- 
tions of the current cosmological parameters, we make 
a global fit to the cosmological observations with the 
publicly available Markov Chain Monte Carlo package 
CosmoMC [25|,|42|]. We assume purely adiabatic initial con- 
ditions, and impose the flatness condition motivated by 
inflation. Our most general parameter space is 

p = (uj b ,uj c ,Qs,T,n s ,n' s ,a s ,A s ), (I) 

where tub = ^b/i 2 and w c = Q, c h 2 are the physical baryon 
and cold dark matter densities relative to the critical den- 
sity, Og characterizes the ratio of the sound horizon to 
the angular diameter distance at decoupling, r is the op- 
tical depth and A s is defined as the amplitude of the 
primordial power spectrum at k = 0.05 Mpc -1 . The 
pivot scales for n s and a s are also chosen at k = 0.05 
Mpc -1 . Here h is the Hubble constant in unit of 100 
km s- 1 Mpc- 1 . The bias factors of 2dF, SDSS and 
LUQAS+CROFT sample of Lyman a (Lya) as analyzed 
in Ref. |40j have been used as nuisance parameters. 

For some of our simulations we have included the fol- 
lowing priors: For h, we make use of the Hubble Space 
Telescope (HST) measurement by multiplying the likeli- 
hood function with a Gaussian centered around h = 0.72 
and with a standard deviation a = 0.08 We impose a 
weak Guassian prior on the baryon density from Big Bang 
nucleosynthesis (BBN) 0: Q, b h 2 = 0.022 ±0.002 (1 a). 
Simultaneously we will also use a cosmic age tophat prior 
as 10 Gyr< t <20 Gyr. 

In our calculations we have taken the total likelihood 
to be the products of the separate likelihoods (A) of 
CMB, LSS, SNIa and Lya. In other words defining 
Xz,,i = -2 log A, we get 

2 _ 2 i_2 |2 i_2 (n\ 

XL,total — Xh.CMB t Xl,LSS ~+~ Xh,SNIa "+" XL,Lya ■ \ z ) 

If the likelihood function is Gaussian, x\ coincides with 
the usual definition of \ 2 up to an additive constant cor- 
responding to the logarithm of the normalization factor 
of C. In the computation of CMB we have included the 
three-year WMAP (WMAP3) data with the routine for 
computing the likelihood supplied by the WMAP team 
|14|. For the purpose of comparisons in some cases we 
also include the first-year temperature and polarization 
data 0, Q with the routine for computing the likelihood 
supplied by the WMAP team [!j . To be conservative but 
more robust, in the fittings to the 3D power spectrum 
of galaxies from the SDSS [2(j we have used the first 14 
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bins only, which are supposed to be well within the lin- 
ear regime |45|. In the calculation of the likelihood from 
SNIa we have marginalized over the nuisance parameter 
|46j. The supernova data we use are the "gold" set of 
157 SNIa published by Riess et al in Ref. [21| and the 
71 high redshift type la supernova discovered during the 
first year of the 5-year Supernova Legacy Survey (SNLS) 
[22j| . respectively. In the fittings to SNLS we have used 
the additional 44 nearby SNIa, as also adopted by the 
SNLS group |23. In order to be conservative but more 
robust, we did not try to combine SNLS together with the 
Riess sample for cosmological parameter constraints 1 . 

The Lyman-a forest corresponds to the Lyman-a 
absorption of photons travelling from distant quasars 
(z ~ 2 — 4) by the neutral hydrogen in the intergalac- 
tic medium (IGM). It is theoretically possible to infer 
the matter power spectrum from the Lyman-a forest 
EE El EE El E3 • Previously Refs. EE || have made 
some studies on the shape of the primordial spectrum and 
<j 8 using Lya data after the release of WMAP3. Here we 
use the LUQAS+CROFT sample of Lyman a as analyzed 
in Ref. E3 ( dubbed V-Lya) and the SDSS Lyman a 
sample from Ref. E2 (dubbed S-Lya) for our study. We 
do not try to combine the two sets of Lya simultaneously 
given the possibly unknown systematics of Lya |55| . 

The likelihood of V-Lya [40] has been included in the 
publicly available package CosmoMC EE E2, which can 
be used directly to probe some scale dependence of the 
primordial spectrum. On the other hand for S-Lya an- 
alyzed in Ref. 41] we need to use the online C++ code 
|56| . The package of CosmoMC being written in Fortran 
90, it is not so straightforward to apply the code in Ref. 
56] in CosmoMC. Recently Slosar [53 made a patch of 
SDSS Lya to CosmoMC. We have corrected the patch E3 
after some nontrivial crosschecks and made our analysis 
for S-Lya basing on the corrected patch. 

Regarding the first-year WMAP data the E mode po- 
larization not being directly measured, there have been 
simultaneously low TT multipoles and high TE multi- 
poles 0, IE 13] • To be conservative Ref. 58] used only 
TT data to probe the possible new physics during in- 
flation and Ref. E3 showed that if assuming without 
significant contaminations, in the general framework the 
observed TE amplitude on the largest scales would be in 
very large discrepancy with the observed TT for the con- 
cordance ACDM model. Moreover Bridle et al. E3] found 
that the claimed preference of a negative a s was merely 



1 We thank Chris Lidman for enlightening discussions on the possi- 
ble systematics led by combining Riess sample L'll directly with 
the SNLS sample |22| in the simple way, and Xiao-Feng Wang 
for stimulating discussions o n th e SNIa correlations as probed by 
the BATM method in Ref. H, MLCS2k2 in Ref. gl, SALT 
in Ref. 49| and th e color parameter AC12 as a new luminosity 
indicator in Refs. loll loll . We also thank Xiao- Feng Wang for 
showing us his preliminary results in combining SNLS with the 
Riess sample in a more professional way |52| . 



due to the lowest WMAP1 multipoles. In order to probe 
the sensitivity of the running to the lower multipoles we 
analyze the running properties using the CMB data with 
and without the contributions of lower multipoles which 
suffer from large cosmic variance |61|. 

For WMAP3 the large scale CMB data not only in- 
clude the temperature power spectrum TT, but also the 
measurements of the polarization spectrum EE, BB and 
the temperature polarization spectrum TE. Given the 
distinctive feature of low I CMB El El> for WMAP3 
the pixel-based method has been applied to Cj T for 
2 < i < 12 and to polarization for 2 < I < 23 in 

Refs. C3 m e m. These large-scale components 

suffer from larger cosmic variance, possibly uncleaned 
foreground contaminations and their statistical nature 
is somewhat different from Gaussian. For these reasons 
we have tentatively prepared two different data-sets of 
WMAP3, one being the "normal" one with the full range 
TT, TE and other polarization data, the other being the 
truncated one where all of the components with I < 23 
have been dropped. We analyze these two data-sets sep- 
arately in various models and compare with each other. 
Also for a clear comparison and for the fact that WMAP1 
should exhibit similar properties on low I to WMAP3, 
we also truncate the TT and TE likelihood of WMAP1 
at I = 24. Thus for the truncated data-set essentially 
for both WMAP1 and WMAP3 we will use only the TT 
and TE CMB spectra on the truncated smaller scales. 
The truncated data-set, although with less amount of in- 
formation, may be more conservative but more robust 
compared with the normal one. 

While the simplest single field inflation generically pre- 
dicts nearly scale-invariant primordial spectrum on the 
CMB relevant large scales, the inflaton would run faster 
in its later epoch and for the scales probed by WMAP+ 
Lya, an appreciable deviation from scale invariance is 
possibly present even in the framework of simple single- 
field inflation models. On the other hand, even for the 
simplest inflaton potentials known for a long time, the re- 
sulting primordial spectrum can be certainly more com- 
plicated than usually expected [3 EE HE E2 • 

Motivated by these we mainly focus on our fittings to 
the Power Law ACDM (PLCDM) Model with a constant 
primordial spectral index n s and Running Spectral Index 
ACDM (RLCDM) Model with a constant running a. In 
order to see if models with more drastic scale-dependence 
are observationally more plausible, we also consider a 
Running Spectral Index ACDM Model with a step-like 
n s divided at k = 0.1 Mpc -1 (SRLCDM) as an example. 
In this model, we assume there is a constant primordial 
spectral index n' s with no running for k > 0.1 Mpc -1 . 
While the amplitudes of n s and n' s are not continuous, we 
naturally assure the continuity of the primordial scalar 
spectrum P s (k) at 0.1 Mpc -1 . 

The parameter erg, being defined as the amplitude of 
density fluctuations at 8 h~ 1 Mpc, cannot be determined 
by WMAP3 alone, which can probe the the shape of the 
primordial spectrum only up to k ~ 0.1 Mpc -1 |10| . with- 
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out extrapolating the information of the primordial spec- 
trum to much smaller scales, except for an exotic case 
that the Hubble constant is anomalously small, which is 
certainly impossible at the first glance in the concordance 
cosmology. On the other hand if we assume NO prior on 
the shape of the primordial spectrum and reconstruct the 
shape of the primordial spectrum in a model indepen- 
dent way like Cosmic Inversion |6.'l (i 1. I65l 66J, l67| , even 
the extremely exotic values of h still remain possible, at 
least in cases we use only the CMB TT observations. In 
the present work we do not pursue the shape of the pri- 
mordial spectrum in a model independent way and will 
study the cosmological implications on the shape of the 
primordial spectrum and the amplitude of as with the 
forementioned method. 

It is believed that the secondary effects on CMB like 
the Sunyaev-Zel'dovich (SZ) effects and CMB lensing ef- 
fects are typically very small and their effects are compa- 
rable with each other. One would take the risk of getting 
biased if considering only one of them [6£|. Neglecting 
both of them is a plausible way in the determination of 
the cosmological parameters for the interest of the cur- 
rent paper 68]. Hence in the present study we have ne- 
glected both of the secondary effects. 

For each regular calculation of MCMC, we run 6 inde- 
pendent chains comprising of 150,000-300,000 chain ele- 
ments and spend thousands of CPU hours to calculate on 
a supercomputer. The average acceptance rate is about 
40%. To get the converged results, we test the conver- 
gence of the chains by Gelman and Rubin 69] criteria and 
typically get R — 1 to be less than 0.05, which is more 
conservative than the recommended value R — 1 <0.1 2 . 



III. ANALYSIS OF CURRENT OBSERVATIONS 
WITH AND WITHOUT LOW MULTIPOLE CMB 
CONTRIBUTIONS 

A. WMAP data alone 

We start with an illustrative figure. In Fig. ^ wc 
delineate the CMB power spectra for the best fit RL- 
CDM Model to normal WMAP3 where the full range 
of the WMAP3 observational data has been useed and 
to the case where the I < 23 data of WMAP3 are not 
used. The dots with error bars are binned WMAP3 
TT and TE data 0. HST, BBN and age priors have 
been adopted. We have cub = 0.0205, uj c = 0.116, 
©s = 1.038, r = 0.092, log[10 10 A s ] (0.05Mpc" 1 )= 3.039, 
n^O.OSMpc- 1 )^ 0.856 and a s (0.05Mpc- 1 )= -0.0634 
for the best fit case with normal WMAP3. Corre- 
spondingly for the case without small I contributions 



2 For the case with WMAP3 only fitting to the Running Spec- 
tral Index ACDM Model with step-like n B truncated at k = 0.1 
Mpc - , which will be shown in the later part of this paper, we 
have R - 1 ~ 0.08. 
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FIG. 1: CMB power spectra for the best fit Running Spectral 
Index ACDM Model to normal WMAP3 where full range of 
the data has been used and to the case where the I < 23 data 
of WMAP3 are not used. The dots with error bars are binned 
WMAP3 TT and TE data. HST, BBN and age priors have 
been adopted. 



the best fit values are Ub = 0.0190, lu c = 0.141, ©g = 
1.037, t = 0.0403, log[10 10 A s ] (0.05Mp C - 1 )= 3.005, 
n s (0.05Mpc- 1 )= 0.747 and a s (0.05Mpc- 1 )= -0.150, re- 
spectively. The values of the minimum \l wm be pre- 
sented later in our Table 5. We can find that the main 
differences of the two results lie on scales I < 20. TE is 
relatively not as good a probe as TT, as can also be seen 
from the small-scale WMAP3 data. 

For the simple PLCDM and RLCDM models some 
of the resulting cosmological parameters from WMAP3 
alone are somewhat different from the previously believed 
concordance cosmology, namely the deviation from the 
scale invariant primordial spectrum and a lower as, as 
well as a mild preference of the running a s . For this pur- 
pose we will investigate the possibility how such discrep- 
ancies are affected if we omit low i contributions which 
suffer from relatively large cosmic variance. In this sec- 
tion no priors are adopted on Hq, Qbh 2 and cosmic age. 

First we focus on the case with the WMAP data 
only. In Table 1 we list the median values and la con- 
strains on PLCDM Model with WMAP1 and WMAP3 
with/without I < 23 CMB contributions, shown together 
with the minimum \\ values. "Normal" stands for the 
cases where the full range of the WMAP observational 
data has been adopted. For simplicity the amplitudes 
of the primordial spectrum A s have not been depicted. 
We should point out that the likelihood of "Normal" 
WMAP1 is somewhat problematic as mentioned in the 
previous parts of the present paper. On the other hand 
we are merely using it for comparison. From Table 1 we 
can find the improvement of precision on the determi- 
nation of the cosmological parameters between WMAP1 
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TABLE 1. Median values and la constrains on Power Law ACDM Model with WMAP1 and WMAP3 and with/without 
£ < 23 CMB contributions, shown together with the minimum \\ values. "Normal" stands for the cases where the full range 
of the WMAP observational data has been adopted. NO priors have been adopted. 





WMAP1 


WMAP3 


No Prior 


Normal 


£ < 23 dropped 


Normal 


£ < 23 dropped 


n b h 2 


0246+ 0026 
u.uz<±o_ 0021 


0.0250 ± 0.0025 


0.0222 ± 0.0007 


0.0223 ±0.0010 


fl c h 2 


0.115 ±0.018 


103+ - 023 


0.106 ±0.008 


0.104 ±0.010 


T 


< 0.449(95%) 


< 0.488(95%) 


< 0.139(95%) 


< 0.293(95%) 


n s 


i n9 +o.os 

L - u ^-0.06 


1.04 ±0.08 


0.955 ±0.016 


9 60+ 033 
u.»ou_ 031 




0.949tr« 


qS9 +0132 

U.t>Ot>_Q 12 7 


0.762 ± 0.050 


0.795 ± 0.065 


H 


75.6+ 9 8 * 




73.0 ±3.2 


73.6 ± 4.9 


xl 


1429.0 


1386.0 


11252.6 


5236.0 



and WMAP3. For the parameter t only the 2a upper 
bounds have been delineated. We should point out the 
reason that our WMAP1 normal results differ signifi- 
cantly from that by WMAP team in Ref. Q is mainly 
due to the differentpriors on r: in our case r is almost 
free while in Ref. Q a stringent prior r < 0.3 has been 
adopted. We can find dramatically for WMAP1 the con- 
straint on t changes very little without the presence of 
low i CMB data. This can be easily understood given 
the fact that for WMAP1 the polarization data of EE 
has not been publicly available. On the other hand for 
WMAP3, all of the constraints have been weaker with- 
out the presence of low £ WMAP and especially for the 
cases of r and erg. For WMAP3 the eminent effects of 
omitting £ < 24 components are that n s is now consis- 
tent with the Harrison-ZePdovich-Peebles scale-invariant 
spectrum and that the value of as is now nontrivially 
higher and less stringently constrained. 

Compared with WMAP1, in WMAP3 the deviation 
from scale-invariance is remarkable due to better mea- 
surements of the polarization and also more sensitive 
measurements of the high £ TT data. We should point 
out our result in Table 1 indicates that at least for the 
relevant scales 24 < £ < 1000, WMAP3 is consistent with 
the scale invariant spectrum at close to la. Thus one may 
conclude the inconsistency of the scale-invariant spec- 
trum claimed by WMAP3 analysis is mainly due to the 
low multipole components. Our constraint of WMAP3 
u £ < 23 dropped" on as is in remarkable agreement with 
the "ACDM ± scale invariant fluctuations (n s = 1)" case 
by the WMAP team on the website 0] , with a consider- 
ably larger error bar. On the other hand, different from 
our case, the previously reported "ACDM ± primordial 
spectrum with sharp cutoff and linear prior" or " ACDM 
± primordial spectrum with sharp cutoff and log prior" 
on the website |14| reported a relatively stringent con- 
straint on n s even in the presence of an additional param- 
eter of the location of the sharp cutoff scale. The reason 
may be partially explained by the fact that the effect of 
reionization is mainly imprinted on low multipoles of po- 
larization data, so that r cannot be determined precisely 



if we omit low £ components 3 . As a result the constraint 
on n s is also loosened as seen in Fig. [21 
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FIG. 2: Two dimensional posterior constraints on n 3 -r 
contours for Power Law ACDM Model with WMAP1 and 
WMAP3 and with/without £ < 23 CMB contributions. The 
solid lines stand for 1- and 2 a respectively. NO priors have 
been adopted. 



In this figure we delineate the two dimensional poste- 
rior constraints on n s -T contours for Power Law ACDM 
Model with WMAP1 and WMAP3 and with/without 
£ < 23 CMB contributions, which correspond to the 
specific subspace of the likelihood surface in Table 1. 



3 We have recovered the results "ACDM + prim ordi al spectrum 
with sharp cutoff and log prior" on the website 14} and also in 
another case we relaxed the prior of "cutoff" , namely we assume 
"ACDM + primordial spectrum with step-like n s at k < 0.05 
Mpc -1 and log prior". We find even in this case a nontriv- 
ial deviation from the Harrison-Zel'dovich-Peebles scale-invariant 
spectrum is still present: n s (0.05 Mpc _1 ) = 0.954 ± 0.0163. 
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The solid lines stand for 1- and 2 a respectively. We 
can see more clearly that for WMAP1 the presence of 
low £ CMB contributions are almost negligible but for 
WMAP3 the contours are significantly different, which 
can be easily understood from the more precise observa- 
tions of WMAP3 and the n s — t degeneracy. 
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FIG. 3: One dimensional posterior constraints on the relevant 
cosmological parameters for Power Law ACDM Model with 
WMAP3 data, with (red) / without (black) I < 23 CMB 
contributions. The red lines stand for WMAP3 constraints 
and the black for constraints with £ > 24 WMAP. NO priors 
have been adopted. 

Correspondingly, in Fig. [3] we show the one dimen- 
sional posterior constraints on the relevant cosmological 
parameters for the PLCDM model with WMAP3 data, 
with/without £ < 23 CMB contributions. The red lines 
stand for WMAP3 constraints and the black for con- 
straints with I > 24 WMAP3. For the sake of clarity 
we have not depicted the WMAP1 cases. One can easily 
find the difference brought forth by the low I WMAP3 
temperature and polarization contributions, especially 
for the constraints on r, n s and as- 

B. WMAP3 data combined with LSS and SNIa 

We now turn to the case LSS and SNIa are com- 
bined with CMB. Here we incorporate both 2dF 0] 
and SDSS HJ as for LSS data, while for SNIa due to 
the forementioned reason we combine in one case with 
the Riess "gold" sample l21l only and in the other case 
with the first year SNLS |22| only. First we study the 
PLCDM model. In Table 2 we show the median values 
and la constrains on the PLCDM model with 2dF, SDSS, 
SNIa (Riess/SNLS sample) and WMAP3 (with/without 



I < 23 CMB contributions), shown together with the 
minimum xj, values. We can easily find that with the in- 
creased observational input, low I CMB components lose 
weight for the less relevant parameters like £lbh 2 , Vt c h 2 
and Hq. On the other hand the parameters like r, n s 
and as are nontrivially affected depending on the pres- 
ence of £ < 23 CMB contributions, because low multipole 
components are much affected by r, while the sensitivi- 
ties of n s and especially as on the \ow-£ components are 
merely due to the fact that we are trying to fit the global 
shape of the power spectrum with less degrees of free- 
dom. Compared with WMAP3 only case in Table I the 
error bars are much smaller in cases of u l < 23 dropped" 
in Table 2. Very intriguingly this leads to the fact that 
now the results are again inconsistent with the Harrison- 
Zel'dovich-Peebles scale-invariant spectrum, and in cases 
one drops the £ < 23 CMB contributions n s is even 
redder (smaller than 1) than the "Normal" case. Also 
for the case of WMAP3+LSS+SNIa(Riess), a smaller 
mean central value of as is present without £ < 23 
CMB contributions, which is in the opposite direction 
compared with WMAP only case. In the left and 
right columns of Table 2 the determinations on the 
cosmological parameters are somewhat different for the 
case with WMAP3+LSS+SNIa(Riess) and that with 
WMAP3+LSS+SNIa(SNLS), but strikingly consistent at 
la. We will address this issue in more details in the later 
part of the current paper. 
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FIG. 4: Two dimensional posterior constraints on n 3 -r con- 
tours for Power Law ACDM Model with 2dF, SDSS, SNIa 
(Riess/SNLS sample) and WMAP3 (with/without t < 23 
CMB contributions). The solid lines stand for 1- and 2 a 
respectively. NO priors have been adopted. 

In Fig. 0] and Fig. [S] we delineate the corresponding 
two dimensional posterior constraints on n s -T contours 
and ag-Q m contours. The solid lines stand for 1- and 
2 a respectively. The left panels stand for cases with 
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TABLE 2. Median values and la constrains on Power Law ACDM Model with 2dF, SDSS, SNIa (Riess/SNLS sample) and 
WMAP3 (with/without I < 23 CMB contributions), shown together with the minimum \\ values. "Normal" stands for the 
cases where the full range of the WMAP3 observational data has been adopted. NO priors have been adopted. 



No Prior 


WMAP3±LSS±SNIa(Riess) 
Normal £ < 23 dropped 


WMAP3±LSS±SNIa(SNLS) 
Normal I < 23 dropped 


n b h 2 
fl c h 2 

T 

n s 
Ho 


0.0222 ± 0.0007 0.0218 ± 0.0007 
0.113 ±0.005 0.113 ±0.005 
< 0.122(95%) < 0.167(95%) 
0.951 ± 0.015 0.937±o.oi8 
0.794 ± 0.035 0.788t£ gte 
69.8 ± 1.8 68.9 ±1.8 


0.0223 ±0.0007 0.0219 ± 0.0007 
0.111 ±0.004 0.111 ±0.004 
< 0.128(95%) < 0.191(95%) 
0.953 ±0.016 0.942 ±0.019 
0.786 ± 0.035 0.786lo.o52 
71.0 ± 1.8 70.1 ±1.8 


xl 


11495.2 5476.8 


11428.4 5410.8 
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FIG. 5: Two dimensional posterior constraints on a$-Q m 
contours for Power Law ACDM Model with 2dF, SDSS, SNIa 
(Riess/SNLS sample) and WMAP3 (with/without £ < 23 
CMB contributions). The solid lines stand for 1- and 2 a 
respectively. NO priors have been adopted. 



WMAP3±LSS±SNIa(Riess) with/without I < 23 CMB 
contributions. The right panels show the correspond- 
ing cases with WMAP3±LSS±SNIa(SNLS). We can find 
that the contours in left panels are consistent but some- 
what different from the right panels, and that a deviation 
from scale invariance is preferred. The constraints on ag 
are less stringent and more consistent, however not as 
significant as the WMAP3 only case as shown in Table 
1. 

For the cosmological implications of 
WMAP3±LSS±SNIa now we turn to the RLCDM 
model. In Table 3 we delineate the median values 
and Id constrains on Running Spectral Index ACDM 
Model with 2dF, SDSS, SNIa (Riess/SNLS sample) and 
WMAP3 (with/without I < 23 CMB contributions), 
shown together with the minimum Xl values. Com- 



pared with Table 2 we can find that almost all of the 
parameters are less stringently constrained with the 
inclusion of a s , especially on the relevant ones like r, 
n s , and as- Due to the inclusion of a s now all of the 
behaviors on the cosmological implications are again in 
the same direction as the WMAP3 only case in Table 1, 
that is, almost all of the error bars are smaller without 
the presence of I < 23 WMAP3 contributions, and a 
larger mean center value of as is present when dropping 
small I CMB contributions. On the constraint on a s 
itself we can find from Table 3 the case of "Normal" 
WMAP3±LSS±SNIa(Riess) favors a negative running 
of a s to nearly 2a. For the other three cases listed a 
running is less favored, but still present and in some 
sense the preference of the negative running is fairly 
independent of low I WMAP3 contributions, which will 
be explored in more details in the remaining part of this 
paper. 

In Fig. we show the corresponding two dimen- 
sional posterior constraints on n s -a contours for Run- 
ning Spectral Index ACDM Model with 2dF, SDSS, SNIa 
(Riess/SNLS sample) and WMAP3 (with/without I < 23 
CMB contributions). The solid lines stand for 1- and 2 
a respectively. We can find that in the case of dropping 
small I CMB contributions the discrepancy between the 
left and the right panels are nontrivial and this needs to 
be reconsidered in more details with the accumulations of 
the currently ongoing SNIa projects. Also the Harrison- 
Zel'dovich-Peebles scale-invariant spectrum lies within 
the 2a range for the case of WMAP3±LSS±SNIa(SNLS) 
without small I CMB contributions, which is different 
from the other three cases. It is also intriguing that in 
the case of "All" combination of WMAP team a neg- 
ative running is preferred to more than 2a, namely, 
a s = -0.061 ± 0.023 0. Their case of "All" differs 
from our analysis here regarding the following aspects: 
they have included small scale CMB observations and 
combined the SNIa data of Ricss and SNLS samples si- 
multaneously, made some detailed analysis on SDSS and 
2dF bias factors and they have considered the secondary 
effects of SZ on WMAP3. And it is interesting to imag- 
ine the consequences without small I CMB contributions 
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TABLE 3. Median values and la constrains on Running Spectral Index ACDM Model with 2dF, SDSS, SNIa (Riess/SNLS 
sample) and WMAP3 (with/without £ < 23 CMB contributions), shown together with the minimum xi values. "Normal" 
stands for the cases where the full range of the WMAP3 observational data has been adopted. NO priors have been adopted 



No Prior 


WMAP3+LSS+SNIa(Riess) 
Normal £ < 23 dropped 


WMAP3±LSS±SNIa(SNLS) 
Normal £ < 23 dropped 


n b h 2 
n c h 2 

r 

n s 

as 
H 


0.0210 ± 0.0009 0.0208 ± 0.0010 

0.116 ±0.005 0.118 ±0.005 

< 0.148(95%) < 0.248(95%) 

0.880l£X 0.862 ± 0.049 

-0 051 2+ 0251 -0 071 7+ 0425 
u.uoi/_ 0254 u.uiii_ 0419 

fl 7OO+0.036 n fin r+0.068 

67.2 ±2.1 66.5 ±2.3 


0.0213 ±0.0009 0.0213 ±0.0010 
0.113±g;gg| 0.114 ±0.005 
< 0.150(95%) < 0.263(95%) 
0.896 ±0.038 0.896 ±0.049 
-0.0422±o.o|^ -0.0484 ± 0.0441 
0.782 ±0.035 0.807i°;^g 
68.9+2' J 68.8 ± 2.3 


Xl 


11491.4 5475.2 


11425.8 5410.2 
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FIG. 6: Two dimensional posterior constraints on n s - 
a contours for Running Spectral Index ACDM Model 
with 2dF, SDSS, SNIa (Riess/SNLS sample) and WMAP3 
(with/without £ < 23 CMB contributions). The solid lines 
stand for 1- and 2 a respectively. NO priors have been 
adopted. 



still significant and the discrepancies on the determi- 
nation of us are rather large among WMAP3 alone, 
WMAP3±V-Lya and WMAP3±S-Lya combinations. In 
cases of dropping small £ WMAP3 contributions, n s is 
nontrivially consistent with scale-invariance. 
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for their case of "All" . 

Although one often takes the risks of uncontrolled sys- 
tematics [55j, Lya can probe scales as small as k ~ lh 
Mpc ~ 1 to yield the amplitude of density fluctuations 
as directly unlike the other data sets. Hence it is urged 
to probe the possible discrepancy between WMAP3 and 
Lya in the determination of the cosmological parameters, 
especially on the value of as- In Table 4 we delineate the 
median values and la constrains on Power Law ACDM 
Model with WMAP3 (with/ without I < 23 CMB contri- 
butions) combined with V-Lya/ S-Lya, shown together 
with the minimum Xl values. One can find again that 
compared with the "Normal" cases of WMAP3±Lya 
combinations where a deviation from scale-invariance is 



FIG. 7: Two dimensional posterior constraints on n s -a$ 
contours for Power Law ACDM Model with WMAP3 only 
(with/without £ < 23 CMB contributions) or WMAP3 com- 
bined with LUQAS±CROFT sample of Lyman a as analyzed 
in Ref. |40j (dubbed V-Lya) or with the SDSS Lyman a sam- 
ple from Ref. | 4l| (dubbed S-Lya). The solid lines stand for 
1- and 2 a respectively. No priors have been adopted. 



In Fig. \7\we show the corresponding two dimensional 
posterior constraints on n s -as contours for the PLCDM 
model with WMAP3 (with/without I < 23 CMB contri- 
butions) combined with V-Lya/ S-Lya. The solid lines 
stand for 1- and 2 a respectively. In the left panel we 
delineate the "Normal" cases and in the right panel we 
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TABLE 4. Median values and la constrains on Power Law ACDM Model with WMAP3 (with/without I < 23 CMB 
contributions) combined with LUQAS+CROFT sample of Lyman a as analyzed in Ref. |40l| (dubbed V-Lya) or the SDSS 
Lyman a sample from Ref. |4j| (dubbed S-Lya), shown together with the minimum \\ values. "Normal" stands for the cases 
where the full range of the WMAP3 observational data has been adopted. NO priors have been adopted. 





WMAP3 + V-Lya 


WMAP3 + S-Lya 


No Prior 


Normal 


t < 23 dropped 


Normal 


I < 23 dropped 


n b h 2 


0.0223 ± 0.0007 


n n99^+ - 0009 
u - uzzo -o.ooio 


0.0228 ± 0.0007 0.0230 ± 0.0010 


Q c h 2 


0.110 ±0.008 


0.109 ± 0.009 


0.120 ±0.006 


0.112 ±0.010 


r 


< 0.139(95%) 


< 0.283(95%) 


< 0.149(95%) 


< 0.313(95%) 


n s 


0.955 ±0.017 


0.958 ± 0.029 


0.964 ±0.016 


n 077+0.032 


as 


0.787 ± 0.047 


0-826™ 


0.857 ± 0.024 


0.891 ± 0.035 


H 


71.6±1;° 


71.8 ±4.4 


68.8 ±2.5 


7i.9±S:° 


xl 


11280.4 


5263.2 


11444.6 


5425.2 



show the corresponding cases where the CMB data arc 
not used for £ < 23. Our left panel is almost exactly 
the same as depicted previously in Ref. [53j], where the 
discrepancies on as are noteworthy for the listed three 
kinds of different data combinations. On the other hand 
as shown in our right panel, in cases we neglect the low 
£ WMAP3 contributions, all of the contours get enlarged 
and these three cases are consistent with each other in a 
wider range of parameters. 



IV. PROBING THE SCALE DEPENDENCE OF 
THE PRIMORDIAL SPECTRUM 

In the above investigations on the discrepancies be- 
tween the implications of the current observations with 
and without low £ WMAP3 contributions we have 
adopted no priors on Hq, £li,h 2 and cosmic age. In the 
following study we mainly focus on the cosmological im- 
plications on the scale dependence of the primordial spec- 
trum, where we prefer to include the physical priors from 
HST, BBN and cosmic age. A running of the primordial 
spectrum will be studied extensively for our remaining 
studies. 



A. WMAP data alone 

We also start from the WMAP only case. In Table 
5 we list the median values and la constrains on Run- 
ning Spectral Index ACDM Model with WMAP1 and 
WMAP3 and with/without £ < 23 CMB contributions, 
shown together with the minimum \l values. As shown 
explicitly for WMAP1 while a negative running of a s is 
close to 2a for the "Normal" case, no running is in full 
consistency with the case u £ < 23 dropped" . Our result is 
consistent with the analysis by Bridle et al. [60|. On the 
other hand, given the problematic likelihood of WMAP 1 
on low £ components, our analysis with "£ < 23 dropped" 
should be more conservative but also more robust com- 
pared with the analysis in Ref. |(3(j ■ With the accumula- 



tion of the observations and better understanding on the 
systematics, the more precise WMAP3 shows a very dif- 
ferent behavior on the shape of the primordial spectrum, 
that is, a negative running is preferred to nearly 2a for 
both of the cases with/without small £ CMB contribu- 
tions. Also the parameter space for as is more consistent 
with the previously well-accepted values for the analysis 
of the cosmic structure formations. 




0.6 0.8 n 1 1.2 0.6 0.8 n 1 1.2 

s s 



FIG. 8: Two dimensional posterior constraints on n s -a s 
contours for Running Spectral Index ACDM Model with 
WMAP1 and WMAP3 and with/without £ < 23 CMB con- 
tributions. The solid lines stand for 1- and 2 a respectively. 
HST, BBN and age priors have been adopted. 

Correspondingly in Fig. |S| we show the resulting two 
dimensional posterior constraints on n s -a s contours for 
the RLCDM model with WMAP1 and WMAP3 and 
with/without £ < 23 CMB contributions. The solid lines 
stand for 1- and 2 a respectively. While for WMAP1 even 
a large positive running is consistent when neglecting the 
small £ CMB contributions, for WMAP3 we can easily 
find the significantly different likelihood space where a 
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TABLE 5. Median values and la constrains on Running Spectral Index ACDM Model with WMAP1 and WMAP3 and 
with/without £ < 23 CMB contributions, shown together with the minimum \\ values. "Normal" stands for the cases where 
the full range of the WMAP observational data has been adopted. HST, BBN and age priors have been adopted. 





WMAP1 


WMAP3 




Normal 


I < 23 dropped 


Normal 


I < 23 dropped 


Q. b h 2 


n n9/is+ ' 0028 
U.U248_ 0027 


0244+ - 0029 
u.uzy^Q 0031 


0.0209 ± 0.0010 


o.oi98i8- 8^ 


n c h 2 


0.0986i°;°l?l 


0.106 ±0.025 


n iic+0.010 

u.iio_ 011 


0.132 ±0.019 


r 


< 0.564(95%) 


< 0.529(95%) 


< 0.155(95%) 


< 0.274(95%) 


n s 




0.996^g 12 7 


o.878±8:82g 


n 78Q+0 088 
u. 1 o»_ 090 


Cts 


-0.08851S 


-0 0432+ 0891 
u.wioz_ 0g86 


-o.o55oiH 3 3? 7 


-0.135 ±0.071 




0.936i° ;£ 


u.aoa_ 134 


n 700+0. 048 


0.831 ±0.063 


H 


82.9tg;< 


78.6±\ 3 3 7 7 


67.7±1;| 


61.4 ±7.0 


xl 


1425.4 


1386.2 


11249.6 


5232.0 



large negative running takes the main parameter space. 
And also for "Normal" WMAP 3 a negative running is 
nontrivially preferred with a higher precision compared 
with WMAP1, where the actual preference should be di- 
minished if the correct likelihood functions is adopted 




B. Combination with other datasets 

For a next step we consider the implications of 
WMAP3 in combination with LSS in great details. In 
Table 6 we show the median values and la constrains 
on the RLCDM model with 2dF, SDSS and WMAP3 
(with/without £ < 23 CMB contributions), shown to- 
gether with the minimum xl values. The separate com- 
binations of WMAP with SDSS and 2dF have also been 
included for the detailed investigations. The discrep- 
ancy between SDSS and 2dF is most significantly de- 
picted in the measurement of the matter density il c h 2 , 
where the tension is even more eminent in cases with- 
out low £ WMAP3 contributions. Our analysis here is 
consistent with the results by the WMAP group [l4j| . 
While WMAP3's better measurements of the TT third 
peak help a lot on the determination of the matter den- 
sity, in the concordance cosmology the measurements 
of fl c h 2 from LSS are also nontrivial. The next re- 
lease of SDSS power spectrum is expected to appear 
soon with much higher precision 4 , the discrepancy on the 
matter density is hopefully to be verified or eliminated. 
Given the current LSS data we find that for the com- 
bination of WMAP3±SDSS the preference of running is 
larger than 2a regardless of the presence of low £ CMB. 
While for the WMAP3±2dF combination or that with 
WMAP3±2dF±SDSS, preference of a negative running 
is less prominent, which nevertheless does exist. 

In Fig. Eland Fig. E3 w e depict the corresponding two 



4 Max Tegmark, private communications. 
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FIG. 9: Two dimensional posterior constraints on n s -a con- 
tours for Running Spectral Index ACDM Model with 2dF, 
SDSS and WMAP3 (with/without I < 23 CMB contribu- 
tions). The solid lines stand for 1- and 2 a respectively. HST, 
BBN and age priors have been adopted. 



dimensional posterior constraints on the n s -a contours 
and on the U m -as contours, respectively. The solid lines 
stand for 1- and 2 a respectively. The left panels show 
the "Normal" constraints and the right ones show the 
cases where low £ WMAP3 are not used for the global 
fittings. 

To measure directly the quantity of erg one needs 
to include the Lya observations. Now we probe the 
scale dependence of the primordial spectrum in the 
light of WMAP3 and WMAP3 in combinations with the 
LUQAS±CROFT sample of Lyman a as analyzed in Ref. 
[iof (V-Lya) or the SDSS Lyman a sample from Ref. 
(S-Lya). We investigate the three cases: PLCDM, RL- 
CDM and SRLCDM cosmology where for SRLCDM we 
mean the Running Spectral Index ACDM Model with 
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TABLE 6. Median values and la constrains on Running Spectral Index ACDM Model with 2dF, SDSS and WMAP3 
(with/without £ < 23 CMB contributions), shown together with the minimum xi values. "Normal" stands for the cases 
where the full range of the WMAP3 observational data has been adopted. HST, BBN and age priors have been adopted. 





WMAP3 + SDSS 
Normal I < 23 dropped 


WMAP3 + 2dF 
Normal £ < 23 dropped 


WMAP3 + SDSS + 2dF 
Normal £ < 23 dropped 


n b h 2 
tt c h 2 
r 

n s 

Ot s 

as 


2 07+ 0010 0194+ - 0012 

u - uzu ' -0.0009 u ' ulyy: -0.0011 

0.123 ±0.008 0.137 ±0.012 
< 0.147(95%) < 0.226(95%) 

0.860+HS 0.767t° :S£ 
-0.0655t° S -0.148«;^ 
0.812 ± 0.040 0.833+° o6o 

64.4 ±3.2 59.114 3 


0.0213 ±0.0009 0.0214 ±0.0010 
0.111 ±0.006 0.112 ±0.007 
< 0.153(95%) < 0.286(95%) 
0.899 ± 0.038 0.8921q.o52 

-0 041 n+ 0256 - 570+ 0467 

U.U41U_ 0258 u - uo ' u -0.0477 

0.7721°;^ 0.812 ± 0.073 
69.8 ±2.5 69.6 ±3.0 


0.0210 ±0.0010 0.0205 ±0.0012 
0.116 ±0.006 0.121 ±0.009 
< 0.149(95%) < 0.254(95%) 

0.882+°;°« 0.845lg;g» 
-0.0505± ;° 2 « -0.0855™ 
0.788 ± 0.036 0.808lo 06° 
67.5l?.; 7 ; 65.4l^ 


xi 


11267.0 5247.6 


11288.8 5273.8 


11307.8 5291.8 



TABLE 7. Median values and la constrains on Power Law, Running Spectral Index ACDM Model and Running Spectral 
Index ACDM Model with step-like n s truncated at k = 0.1 Mpc _1 (SRLCDM), as measured by WMAP3 in combination with 
LUQAS±CROFT sample of Lyman a as analyzed in Ref. Illl (dubbed V-Lya) or the SDSS Lyman a sample from Ref. |4l| 
(dubbed S-Lya), shown together with the minimum xi values. HST, BBN and age priors have been adopted. 







PLCDM 




RLCDM 




SRLCDM 






WMAP3 only 


V-Lya 


S-Lya 


V-Lya? 


S-Lya 


WMAP3 only 


V-Lya 


S-Lya 


10 2 Q b h 2 


2.22 ±0.07 


2.24 ±0.07 


2.27 ±0.06 


2.19 ±0.07 


2.25 ±0.06 


2.12 ±0.09 


2.12 ±0.09 


2.12 ±0.09 


Q c h 2 
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0.110 ±0.007 
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< 0.136(95%) < 0.138(95%) < 0.147(95%) 


< 0.152(95%) 


< 0.157(95%) 


< 0.157(95%) 


< 0.162(95%) 


< 0.170(95%) 


n s 


0.953 ±0.016 


0.956 ±0.015 


0.964 ± 0.015 


0.923i°;° 2 2 6 7 


0.946 ± 0.020 


0.887 ± 0.042 


n QQO+0 039 
U.88^_ .040 


+0.039 
U.O( O_o.038 


n' s 












No constraint 


0.929 ± 0.090 


0.958 ± 0.063 
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Ho 
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preferred at > 2a on scales k < 0.1 Mpc -1 for both of 
the cases and the implications on as turn out to be more 
consistent. This should be relevant to the fact that the 
introduction of Lya data changes the preference of Q m 
compared with WMAP3 only case, where fi m is much 
relevant to the TT third peak. A more intriguing as- 
pect comes from the determination on n' s for the cases in 
combinations with V-Lya or S-Lya: for both cases the 
Harrison-Zel'dovich-Peebles scale- invariant spectrum fits 
well the observations on scales k > 0.1 Mpc -1 . On the 
other hand compared with a constant running RLCDM 
case, a constant negative running for k < 0.1 Mpc -1 
and nearly scale invariant spectrum for scales k > 0.1 
Mpc" 1 is preferred at close to 2a for the WMAP3±S-Lya 
combination and similar behavior also exhibits for the 
WMAP3±V-Lya combination, which is relatively weaker 
than the former case. 

In Fig. Illl we depict the two dimensional posterior 
constraints on Q, m -as contours for Power Law, Running 
Spectral Index ACDM Model and Running Spectral In- 
dex ACDM Model with step-like n s truncated at k = 0.1 
Mpc -1 , as measured by WMAP3 alone and WMAP3 in 



step-like n s truncated at k = 0.1 Mpc -1 , with n' s being 
the additionally introduced parameter which represents 
the spectral index on scales k > 0.1 Mpc -1 . Note in 
this case we have assumed a constant n' s . In Table 7 
we list the median values and la constrains on the de- 
terminations of the relevant cosmological parameters for 
the three cases. For the PLCDM case a deviation from 
scale-invariance is present for all of the three cases in the 
left column of Table 7. For the RLCDM case as we are 
assuming a constant running on all scales, a vanishing 
running is consistent in either combination of WMAP3 
with V-Lya or S-Lya. For the case of SRLCDM we find 
many interesting implications. Because WMAP3 cannot 
probe scales with k > 0.1 Mpc -1 , we cannot get any con- 
straint on the parameter n' s from WMAP3 alone. And 
for this case we get a% < 2.36 at 95% confidence level, 
which we will explain later in more detail. Comparing 
Table 7 with Table 5, we can find that the SRLCDM 
case of WMAP3 only is similar to that in Table 5 for 
the RLCDM case only, and a running is favored at less 
than 2a. Interestingly for the cases of WMAP3 in com- 
binations with V-Lya or S-Lya, a negative running is 
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FIG. 10: Two dimensional posterior constraints on O m -crg 
contours for Running Spectral Index ACDM Model with 2dF, 
SDSS and WMAP3 (with/without I < 23 CMB contribu- 
tions). The solid lines stand for 1- and 2 a respectively. HST, 
BBN and age priors have been adopted. 



combination with the LUQAS+CROFT sample of Ly- 
man a as analyzed in Ref. [40j (dubbed V-Lya) or the 
SDSS Lyman a sample from Ref. j^H (dubbed S-Lya). 
We can easily find that the discrepancy is the largest 
in PLCDM and as we introduce the additional param- 
eters like a constant running and also then work in the 
SRLCDM framework, they gradually become consistent 
with each other. And for the case of SRLCDM, as n' s 
is not constrained, the resulting constraint on eg would 
definitely be extremely weak and we put the WMAP3 
constraint alone on the right panel. 

For a better study on the cosmological constraint by 
WMAP3 alone, in Fig. ^] we show the one dimensional 
posterior constraints on the relevant cosmological param- 
eters for the PLCDM model, the RLCDM model and the 
SRLCDM model respectively, as measured by WMAP3 
alone. The black lines stand for PLCDM, red lines for 
RLCDM and the blue for SRLCDM model. We should 
point out that although our chains have converged, our 
results on the SRLCDM model is much dependent on 
the prior adopted on the parameter n' s , which is easily 
understood given the fact that n' s is unconstrained for 
WMAP3 alone. And a weaker prior on n' s would lead to 
a weaker constraint on ag- 
in Fig. EH we delineate the resulting 2a posterior con- 
straints on n s (k) for Running Spectral Index ACDM 
Model and Running Spectral Index ACDM Model with 
step-like n s truncated at k = 0.1 Mpc , measured by 
WMAP3 in combination with LUQAS+CROFT sample 
of Lyman a as analyzed in Ref. |40| and the SDSS Ly- 
man a sample from Ref. |4l| . One can easily find the 
scale dependence on either of the upper or the bottom 



figures for the case of SRLCDM model. 



V. DISCUSSION AND CONCLUSION 

Through our detailed investigations we have found 
that values of cosmological parameters obtained by 
MCMC change considerably depending on whether we 
use low multipole components of CMB or not. Cur- 
rently, the WMAP team use a very large number of pixels 
(666+1172 pixels for low multipole temperature and po- 
larization spectra 0, HTl Il2l 0|), which may take too 
large a weight for the determination of the likelihood and 
hence cosmological parameters. This may be somewhat 
problematic because these low multipoles suffer from rel- 
atively large cosmic variance. On the other hand, if we 
incorporate low £ polarization components, the optical 
depth r can be determined with a smaller ambiguity be- 
cause the degeneracy between n s and r is partially re- 
moved, which in turn result in a better determination of 
the spectral index. Hence it is difficult to conclude which 
approach gives more precise results. 

We obtain different values of cosmological parameters 
with different datasets. For example, we find a very low 
H in the combination of WMAP3 (without I < 23) and 
SDSS is rather interesting, even though we have adopted 
the HST prior in Table 6. Nonetheless these discrepancies 
are typically within 2a and do not imply the datasets we 
have used are incompatible with each other, including the 
differences between the Riess "gold" sample and the first 
year SNLS, SDSS and 2dF, or LUQAS+CROFT sample 
of Lyman a as analyzed in Ref. [40l | and the SDSS Lyman 
a sample from Ref. 41] 

While our results as a whole support the concordance 
cosmology, they also indicate we are still far from being 
able to determine all of the six cosmological parameters 
with two digits' accuracy in the context of this simple 
power-law ACDM model. Hence we should be open- 
minded to new physics which may be hidden in the error 
bars of the current observational data and continue blind 
analysis without any theoretical prejudices. A good ex- 
ample is the possible running of the spectral index, which 
is preferred to more than 2a in a number of combinations 
of datasets. 
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